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Abstract 

We compute the shear viscosity of a two dimensional Fermi gas interacting via a short range 
potential with scattering length a2d in kinetic theory. We find that classical kinetic theory predicts 
that the shear viscosity to entropy density ratio of a strongly interacting two dimensional gas is 
comparable to that of the three dimensional unitary gas. We apply our results to the damping of 
collective modes of a trapped Fermi gas, and compare to experimental data recently obtained in 
E. Vogt et al., larXiv:llll.lT73l 
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I. INTRODUCTION 



The study of transport properties of strongly interacting, scale invariant or approximately 
scale invariant fluids has led to many recent discoveries that connect the physics of cold 
atomic gases, properties of the quark gluon plasma, and quantum gravity jl]. Nearly ideal 
hydro dynamic flow in cold atomic gases was observed in the expansion of a dilute Fermi 
at unitarity [2j, and similar results were observed in heavy ion collisions at the relativistic 
heavy ion collider (RHIC) [3J . Recent analyses show that both the quark gluon plasma and 
the dilute Fern. gas at unitarity a r e characterised by a shear viscosity to entropy density 

|. This result is close to the value r]/s = h/ (Anks) which was found 



ratio r]/s < 0.5 h/k B (4 
in the strong coupling limit o 



a large class of field theories that can be analyzed using the 



AdS/CFT correspondence [§[ ^]. 

The AdS / CFT result is independent of the dimensionality of the fluid, and it is interesting 
to study whether nearly perfect fluidity can be observed in two-dimensional fluids. It was 



suggested, for example, that electrons in graphene might behave as a nearly perfect fluid [10 ]. 
Recently, a group at the Cavendish Laboratory investigated the damping of collective modes 
in a cold atomic Fermi gas tightly confined in one direction ll|. Vogt et al. determined the 
damping constant as a function of T/Tp in the range T /Tp = (0.3 — 0.8), and for different 
interaction strengths \og(kpa 2 d) = (2.7 — 42). Here, T/Tp is the temperature in units of the 
Fermi temperature, kp is the Fermi momentum, and a 2 d is the two-dimensional scattering 
length. In the present work we compare these results with the predictions of kinetic theory. 
Formally, kinetic theory is reliable in the limit of high temperature, T ^> Tp, or in the case 
of weak interactions, Kpa 2 d ^> 1. In the case of the three dimensional Fermi gas at unitarity 

icability of kinetic theory is larger than one might 
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it was observed that the range of app 
expect, extending down to T ~ 0.4 Tp 

II. KINETIC THEORY 

The viscous stress tensor in hydrodynamics is given by SHy = —rjcTij — (Sij{a) with 

2 

o-y = diVj + djVi - - 5ijd k v k , (1) 

and (ct) = dtv Here, Vi is the flow velocity and d = 2,3,... is the number of spatial 
dimensions. We will determine 7] by matching the hydrodynamic result to kinetic theory. 



The stress tensor in kinetic theory is given by 



(2) 



where v is the number of degrees of freedom (v = 2 for a two-component Fermi gas), dT p = 
is the volume element in momentum space, and 5f p is the off-equilibrium correction to 



the distribution function. We will use the ansatz 



Sf P = fp - jrXij(p)<7i. 



\ </(/') =PijX 



( P 2 \ 

;Y - — = 



(3) 



\2mT) 

where f p is the classical equilibrium distribution function and p^ = piPj — \bijP 2 - We will 
study the role of quantum statistics below. We compute Sf p by solving the Boltzmann equa- 

2 

tion for fermions with dispersion relation E p = subject to elastic two-body scattering. 
At this level of approximation the bulk viscosity vanishes. This is the correct result for 



3d fermions at unitarity 
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-Il7j. but the bulk viscosity is expected to be non-zero for 3d 



fermions away from unitarity, and for 2d fermions at any value of the scattering 



18- 



ength. The 



20|, but the 



dependence of the integral / du)((u}) on [kpo)~ l is constrained by sum rules 
bulk viscosity at zero frequency has not been determined. Vogt et al. measured the damping 
of a 2-d quadrupole mode [U], which is not sensitive to bulk viscosity. 

Matching the kinetic theory expression for 5Uij to hydrodynamics we get 

2u 1 



V 



{Pij I Xij ) 



(4) 



(d-l)(d + 2) mT 

where we have defined the inner product (a\b) = J dT p f p a(p)b(p). The function Xij(p) is 
determined by the linearized Boltzmann equation 



7^ M = C\Xij) ■ 



(5) 



Here C is the linearized collision operator C\xij) = \C[Xij}) with 

C\xM\ = n(/ dT *) f°(P2)(2n) d+1 5 d (P - P')8{E - E')\T\ 2 

■ [Xiji.Pi) + Xij(p2) - Xijips) ~ Xij(P4>] 



(6) 



where T is the T-matrix for elastic two-body scattering 12 — > 34. We have also defined 
Pi j2 = y ± q, ps,4 = -y ± q', E = E P1 + E P2 and E' = E P3 + E P4 . Given the T-matrix we can 
determine Xij from equ. (jSJ) and then compute the shear viscosity using equ. (jl]). In practice 
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it is useful to reformulate the calculation as a variational problem. The shear viscosity can 
be written as 



1 



V 



(XijlPij)' 



(7) 

follows from 



(d-l)(d + 2) m 2 T (xy|C|x«> ' 
The equivalence of this result and the previous expression given in equ. 
the linearized Boltzmann equation. The result is variational in the sense that for a trial 
function x™ r equ. (J7J) provides a lower bound on the shear viscosity. The exact solution 
of the linearized Boltzmann equation can be found by maximizing equ. (jTJ). In the three 
dimensional case it is known that the quadratic ansatz Xij = Ma is an excellent solution, 



2 If ] . We will see that despite 



providing results for the shear viscosity that are accurate to 2% 
the different structure of the scattering amplitudes in two and three dimensions the matrix 
elements of the collision operator are very similar. We will therefore use the trial function 

Xij Pij- 

In two dimensions the scattering matrix for elastic scattering mediated by a short range 
potential is given by [22] 

T = i — 7 9 2 \ i • ' (8) 

m - logtya^j + itc 

where a 2 d is the two-dimensional scattering length. The cross section is = ^-|T| 2 . The 
matrix element of the linearized collision operator can be reduced to a one-dimensional 
integral. We find 



(Xii\C\Xii) = 4T(mT) 3 / dx — T —— l - 

V " ' 3 ' 1 ^ log 2 (a; 2 T/T aiM )+7r 2 



(9) 



where we have defined T a 2 d = l/( mo 2d)- The integral in equ. © can be computed using 
the saddle point approximating. This amounts to replacing the term x 2 in the denominator 
by 5/2. The final result for the shear viscosity is 

2 



V2d 



mT 
2^2 



log 



5T 



2T a . 



2d . 



(10) 



where we have set u, the number of spin states, equal to two. We can use this results to 
compute the dimensionless quantities r]/n and rj/s. We find 

l 2N 



V2d 

n 




log 



5T 



IT, 



a,2d. 



(11) 



where Tp c = (kp c ) 2 /(2m) is a function of the local Fermi momentum, kp c = (27m) 1 / 2 . The 
entropy per particle is s/n = log(T/Tj? c ) + 2. 
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FIG. 1: The left panel shows the viscosity to density ratio tj/n as a function of T/Tp for a two- 
dimensional Fermi gas with (kpa) 2 = 2. Here, Tp = k F /(2m) and &p = (27m) 1 / 2 characterize the 
homogeneous Fermi gas. The solid line includes the effects of quantum statistics, the dashed line 
shows the high temperature limit given in equ. (jlip . and the dotted line shows the low temperature 
limit. The right panel displays the shear viscosity to entropy density ratio. The dash-dotted line 
shows the proposed bound rj/s = 



It is instructive to compare these expressions to the analogous formulas in three dimen- 
sions. The T-matrix is 



r 



4tt 1 



m -a 3d + iq 



(12) 



and the cross section is ^ = ^2 17~| 2 . The collision integral is 



are 



16m 7/2 T 9/2 f°° 
(Xv\C\Xij) = 37r 5/2 J dx \ + TaM /(x*T) 



(13) 



where T a ^d = l/( ma :L)- At unitarity, — >■ oo, the integrand differs from the result in two 
dimensions only by logarithmic terms. The shear viscosity at unitarity is 

15 



32^ 



(mT) 



.3/2 



(14) 



In the limit T a ^/T > 1 we find t]^ = 5(mT) 1 / 2 /(32 v /7ra 2 ). The three dimensions the 
density is n = (k l p c ) 3 / (37r 2 ) , and the shear viscosity to density ratio is 



45tt 3 / 2 / T 



V3d 



n 64x72 V t f 



3/2 



loc 



(15) 



Finally, the entropy per particle is s/n = | log(7rT '/Tp) + log(3/4) + 5/2. 
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The results in two dimensions are plotted as the blue dashed lines in Fig. [TJ We have 
chosen {kpa2d) 2 = 2, which means that the two body binding energy Eb = l/(ma| d ) is equal 
to the Fermi energy. This corresponds to the BEC/BCS crossover regime. We observe that 
for T/Tp c < 0.5 the shear viscosity to entropy density ratio reaches f]2d/ s — 0.5, comparable 
to the result for the three dimensional Fermi gas at unitarity. In this regime kinetic theory 
is not reliable - effects due to quantum statistics, correlations and fluctuations are likely 
to play a role. Quantum statistics can be included straightforwardly in the kinetic theory 
calculation by including appropriate statistical factors in equ. (j2|3|) and ()6]). The result is 
shown as the solid line in Fig. [TJ Pauli blocking suppresses the scattering matrix element 
and leads to r] 2 d/n ~ (Tp c IT) 2 [log (Tjf c (T) } 2 as T — > 0. This result is expected from Landau 



Fermi liquid theory 



231 ] . We observe that in two dimension the effect of Pauli blocking 



is quite large, but we also emphasize that at strong coupling the inclusion of quantum 
statistics is not necessarily an improvement over the classical calculation. In the case of 
thermodynamic quantities, like the second Virial coefficient, it is well known that effects of 
quantum statistics appear at the same order in T/Tp as higher order interaction terms. A 
similar effect is seen in the many body T-matrix calculation of the shear viscosity of the three 



dimensional gas at unitarity by Enss et al. 



16|. These authors include pairing correlations 



and vertex corrections in addition to the effects of quantum statistics. They find that the 
shear viscosity to entropy density ratio remains very close to the classical result even in the 
very degenerate regime T ~ (0.2 — 0.5) Tp. 



III. DAMPING OF COLLECTIVE MODES IN A TRAPPED GAS 

In hydrodynamics the damping of collective modes is governed by the rate of energy 
dissipation 

E = -^Jd 3 x V (x)(a lJ ) 2 , (16) 

where we have neglected bulk viscosity and assumed that the system remains isothermal 
(so that heat conductivity can be neglected). For simple modes like the quadrupole oscil- 
lation studied by Vogt et al. the velocity field is linear in the coordinates and the stress 
tensor is spatially constant. In this case the decay rate is sensitive to the spatial integral 
of rj(x). On dimensional ground we can write the viscosity of the homogeneous system as 
r] = na n (T/Tp c , k l p c a). The spatial integral over r](x) can then be written as N(a n ), where 
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N is the total number of particles and (a n ) is the value of a n averaged over the density dis- 
tribution of the cloud. In the hydrodynamic regime measurements of the damping constant 
of collective modes can therefore be interpreted as measurements of {a n ). 

The difficulty with this approach is that in kinetic theory rj = na n is independent of 
the density and the spatial average (a n ) is formally infinite. Physically, this problem is 
related to the fact that for any finite collective mode frequency hydrodynamics cannot be 
applicable in the dilute corona of the cloud, so that the integral in equ. ffl6|) has to be cut off 
at low density [I4J . In kinetic theory this can be done by taking into account the frequency 
dependence of the shear viscosity 

77(0) 



1 + 



(17) 



where tr is the viscous relaxation time, which is the time it takes for the stress tensor to 
relax to the Navier-Stokes form SUij = — 77(0)0^. We will see that the relaxation time is 
inversely proportional to the density, and that the spatial integral over rj(u) is therefore 



finite 
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l4. 



The relaxation time can be determined in various ways, for example by solving the lin- 
earized Boltzmann equation in a time-dependent velocity field 24], |25(, by computing the 
viscosity spectral function 26], or by evaluating the relaxation time in second order hydro- 



dynamics 



27j . The relaxation time is also constrained by viscosity sum rules 18N20|. Using 



the methods described in 



26j we can show that in kinetic theory t](u>) satisfies the sum rule 

If P 
- / dur](u) = — , 

7T J L 



(18) 



where P is the pressure. This sum rule is valid in both two and three dimensions. Combining 
equ. (|T71) with the viscosity sum rule equ. (|T8l) we get tr = rj/P ~ i]/(nT). 

We note that the sum rule in equ. ( fl8l) follows from the definition of the stress tensor in 
kinetic theory, see equ. (T5]). If the stress tensor is defined as an operator in the quantum 
;heory one finds that the spectral function in two dimensions has al/u tail at high frequency 
2p| . The corresponding behavior in three dimension is p(co) ~ Xj^JZj. This tail does not 
appear in kinetic theory because kinetic theory is an effective theory for energies u < T. 
In the quantum mechanical sum rule the high frequency has to be subtracted. In the high 
temperature regime, T > Tp, the conclusion is the same as before: the high frequency tail 
does not contribute to the sum rule, and the width of the transport peak is controlled by 
the relaxation time r R = i]/(nT) [161 ] . 




FIG. 2: This figure shows the trap average of the shear viscosity to density ratio a n as a function 
of T/Tp for different values of logikF^d)- Here, Tp = ui±N 1 / 2 is the Fermi temperature and 
kp = (2mTp) 1//2 is the Fermi momentum in the trap. The scale is set by oo±, the transverse (two 
dimensional) confinement frequency. We have used N = 4 • 10 3 . We compare our results to the 
data from Vogt et al. [111 ]. 



We can now compute the trap average of T]{oj). We will use the high temperature ap- 
proximation for the cloud density. This is consistent with the classical kinetic calculation of 
rj, and is expected to be a good approximation in the regime T/T F > 0.3 studied by Vogt 
et al. In this limit the density profile of a 2-dimensional cloud is 

mT (T F - 2 



n(x) 



2tt V T 



2 2 ' 

mm, x 



cxp 



IT 



(19) 



where Tp = co^N 1 ^ 2 is the Fermi temperature of the trapped gas. For the 2-dimensional 
quadrupole mode the frequency is given by to = y/2co± 28l-l30|. We note that the quadrupole 
mode is volume conserving, and the frequency is independent of the equation of state. We 
get 



(a n ) = —R ( — ) log 
v ' 2tt \T F ] 



« Ntt z , , 



Tp 
T 



R 



log 



5T 



1 2 



2T, 



a,2d. 



(20) 



This result is plotted in Fig. |2l We observe that for small values of log(kpa2d) and T/Tp 
the trap average (a n ) grows approximately as T 2 . This power law can be understood as one 
factor of T arising from the temperature scaling of rj, and one factor of T from the inverse 
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density at the center of the trap. For larger values of \og(kFCi2d) the growth of the relaxation 
time compensates the growth in 77 and the trap average (a n ) is only weakly temperature 
dependent. 

In Fig. [2] we also compare our results to the data obtained by Vogt et al. 11] . We 
observe that the predicted dependence of (a n ) on T/T F and \og(k F a 2 d) is in qualitatively 
agreement with the data. The theoretical predictions are in quantitative agreement with 
the data for log(/cpa 2 d) = 5.3 and 9.7. The disagreement between theory and data for 
\og{k,Fa2d) = 2.7 is somewhat puzzling, because this value of log(fcpa2d) corresponds to a 
more strongly interacting fluid, and we would expect hydrodynamics to work better. Of 
course, the kinetic theory calculation of the shear viscosity might break down at strong 
coupling and T/T F < 1. Another possible issue is that the experimental analysis used a 
free Fermi gas model to estimate the energy of the mode. At strong coupling this approach 
will tend to overestimate the energy, and the extracted trap average (a n ) is too large. The 
theory also under-predicts the data for large values of \og{kF(i2d)- This is less surprising, 
because hydrodynamics is expected to break down in this regime. 



IV. OUTLOOK 



The observed qualitative agreement between experiment and the predictions of kinetic 
theory suggests that the shear viscosity of the two dimensional Fermi gas can be extracted 
from measurements of the damping of collective modes. In order to do this quantitatively 
a number of effects will have to be studied more carefully. We observe, in particular, that 
for \og{k,Fa2d) ^ 5 the measured collective mode frequencies are not close to the hydro- 
dynamic predictions. This implies that dissipative effects are not accurately described by 
the hydrodynamic expression given in equ. ( Fl6l) . A more appropriate approach is to treat 
the collective mode itself in kinetic theory. This calculation will also provide an indication 
whether the observed damping at large \og{k F a2d) is related to collisions, or other effects 
that are not taken into account in a kinetic or hydrodynamic treatment. 

We note that even though the observed trap averaged values of (a n ) are on the order of 1 
or larger the corresponding value of rj/s at the center of the trap could be quite small, on the 
order of rj/s ~ 0.5, see Fig. [TJ In the interesting regime T < 0.5 Tp classical kinetic theory 
is not reliable. In two dimensions, in particular, correlations and fluctuations are likely to 
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play an important role. An important example of a correlation effect is the pseudo-gap 
phenomenon which was argued to play an important role in the transport behavior of the 
three dimensional gas 31] . A pseudo-gap has been observed in the two dimensional gas in 
the regime \og(kpa2d) 1 (22|- The phase transition in two dimensions is of Berezinsky- 
Kosterlitz-Thouless (BKT) type, and the two dimensional Fermi gas may provide a very 
clean system to study transport properties near the BKT transition. It is also known that 
in two dimensions hydrodynamic fluctuations lead to a slow, logarithmic, divergence of the 
shear viscosity with the system size [33 1. 
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